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The dependence of the energies of axially symmetric monopoles of magnetic 
charges 2 and 3, on the Higgs self-interaction coupling constant, is studied numer- 
ically. Comparing the energy per unit topological charge of the charge-2 monopole 
with the energy of the spherically symmetric charge-1 monopole, we confirm that 
there is only a repulsive phase in the interaction energy between like monopoles 
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It is well known that like monopoles of the Georgi-Glashow (GG) model exhibit only a 
repulsive phase. In the Bogomol'nyi-Prasad-Sommerfield (BPS) limit |3|,[|], when the dimensionless 
coupling constant A parametrising the strength of the Higgs self-interaction potential vanishes, this 
interaction disappears. This is a consequence of the vanishing of all components of the stress tensor 
H in this limit. In the BPS limit the Higgs field becomes massless and mediates a long range 
attractive force which cancels the long range repulsive magnetic force of the U(l) field, exactly. The 
force between like monopoles in the BPS limit was studied in detail by Manton || and Nahm 0, 
who showed that it decreases faster than any inverse power. When A > however, the Higgs field 
becomes massive and as a result decays exponentially. Consequently the long range magnetic field 
dominates at large distances, leading to the repulsion of like monopoles of the GG model. This 
was concluded by Goldberg et al. using the time-rate of change of the stress tensor for the field 
configuration of two exponentially localized monopoles, situated apart at a distance much larger 
than the sizes of the monopole cores. 

Inspite of the well known scenario described above, to our knowledge, there has not been any 
detailed numerical verification of it to date. It is the purpose of the present work to supply this. 

A numerical study of the energy of the unit charge spherically symmetric 't Hooft-Polyakov 
monopole was carried out by Bogomol'nyi and Marinov long ago where the dependence of the 
energy on A was studied in detail. No such study was carried out for axially symmetric monopoles 
of higher charges. The study of axially symmetric monopoles was confined to the BPS limit only 
pOpl] ]. In the present work we will analyse the A > case, and in particular for the n — 2 and 
n = 3 monopoles of the GG model. 

Our procedure is analogous to the work of Jacobs and Rebbi where they study the energy 
per unit charge of the topologically stable charge-n vortices of the Abelian Higgs model ( AHM) ^ , 
and demonstrate that in these theories there are bound states of vortices or more generally, that 
there are both attractive and repulsive phases. 

In both the AHM and the GG model at hand, the masses of the solitons, the vortex (resp.) the 
monopole, depend on the strength of the dimensionless Higgs self-interaction coupling constant A. 
In both cases, topological stability is guaranteed if A > 0. The difference between these two systems 
is that the values of A for which the topological inequalities are saturated, namely the critical values 
A c , are different. For the AHM A c = 1 while for the GG model A c = 0. It is clear that in the BPS 
limit (A = A c ) the masses per unit charge of charge-n solitons in cither of these two theories are 
equal to 1, for all n. 

It follows that in the case of the AHM the energy per unit charge vs. A curves for the n = 1 
and n — 2 vortices cross at least at one point, namely at A c = 1, and hence that there will be both 
attractive/repulsive phases according as to whether the mass per unit charge of the n = 2 vortex is 
lower/higher than that of the n = 1 vortex. That there is only one such intersection point can be 
speculated on the grounds that we expect such curves for solitons to be monotonic. This scenario 
was confirmed by the numerical analysis of Ref . fl2]| . 

In the case of the GG model, the n = 1 and n = 2 curves meet at A c = 0. Speculating again that 
that there is only one such intersection, we would conclude that there is only an attractive/repulsive 
phase according as to whether the mass per unit charge of the n = 2 monopole is lower/higher than 
that of the n — 1 monopole. Speculating further, on the basis of the results found in Refs. 1^,0, we 
would expect that there is only a repulsive phase in the GG models. This expected result is precisely 
what we will verify by numerical analysis in the present work. 

In addition to confirming this expected conclusion, we will find the detailed dependences of the 
monopole masses per unit charge of the axially symmetric charge-n monopole on the Higgs coupling 
constant A. The numerical analysis will be carried out for n = 2 and n = 3. In this sense, the 
present work is a charge-n version of the work of Ref. Q . 

The static Hamiltonian of the GG model is 

H = l -Tr{\F l0 1 2 ) + Tr(\ D^) + ^A(2Tr(0 2 ) + n 2 ) 2 , (1) 
with the field strength tensor 
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F i j=a i Aj-d j A i + [A i ,Aj] (2) 

of the gauge field Aj, and the covariant derivative 

Di<t> = di<f>+ [A h ct>] , (3) 

of the Higgs field </>, and r\ denotes the vacuum expection value of the Higgs field. The mass of the 
vector fields is Mw = V (we have set the gauge coupling constant to one), and the Higgs coupling 
constant A is related to the Higgs mass Mh by A = 2M|/M^. 

Employing spherical coordinates, we parametrize the non-vanishing components of the gauge field 
by the Ansatz 



■ (") ■ (n) 

A v = -nsine[H 3 (r,6) ^- + (1 - H A (r,0)) -^-] . (4) 



and the Higgs field by 



■ (f») ■ (n) 

IT IT 

<P = Mr,O)^- + Mr,0) • (5) 
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The matrices are defined by 



T^ n ' = sin 9 cos n<^ ci + sin 9 sin 02 + cos 9 03 



Tg = cos 6* cos rup <j\ + cos sin nip o 2 — sin 9 03 



where cr; denote the Pauli matrices and the integer n refers to the winding number which equals 
the topological charge of the solutions. This Ansatz (based on the Ansatz used in Ref. (h]]) allows 
for axially symmetric solutions with n > 1 and reduces to the Ansatz for spherically symmetric 
solutions if n = 1. Due to a local U(l) invariance of the Ansatz it is necessary to impose a gauge 
fixing condition which guarantees the uniqueness of solutions. Here we take the gauge condition 
rd r H x -d e H 2 = @. 

The energy of a solution with winding number n 



E n (X) = J Hd s r > 4nQr) (6) 

is bounded from below by the magnetic flux 

£i jk J TriFijDk&Sr = e ijh J Tr{$F l0 )dS k = 4ttQ v , (7) 

where Q represents the topological charge, and Q = n for the Ansatz given above, provided that the 
solutions satisfy the appropriate boundary conditions (stated in the following). 

The boundary conditions follow from the requirements of finite energy and analyticity as well as 
symmetry. They are given by 

ffi(O,0) = .ff3(O,0) = 0i(O,0) = 02(O,0)=O, H 2 (0,e) =H A (0,9) = I 
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at the origin, 

H i {oo,9) = cf> 2 (oo,6)=0 , <(> 1 (oo,6)/ri = l 

at infinity, 

fli(r,0) =H 3 (r,0) = fa(r,Q) = , 9 e if 2 (r,0) = d e H 4 (r,0) = fy&(r,0) = 
on the z-axis and 

ffi(r,7r/2) = fl 3 (r,7r/2)=0 2 (r ) 7r/2) = O , d e H 2 (r,7r/2) = d e H A {r,n/2) = ^(r.Tr/2) = 
on the p-axis. 
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Fig. 1 

The normalized energy per unit charge as a function of A is shown for n = 1 (dashed line) and 
n = 2 (solid line). The inlet shows the excess of the energy per unit charge for small values of A for 
n = 1 (*), n = 2 (+) and n = 3 (x). The straight lines with slope 1/2 are fitted to the data values 
corresponding to A = 1 x 1CP 5 for n = 1 and n — 2 and to A = 2 x 10~ 5 for n = 3. 

The partial differential equations for the functions Hi(r, 9) and <j)i(r, 9) are obtained by inserting 
the Ansatz into the Euler-Lagrange equations of the static Hamiltonian with gauge fixing term 
added. We have solved the partial differential equations numerically, subject to the above boundary 
conditions for the winding numbers n = 1 — 3 and for a large range of values of the Higgs coupling 
constant A. 

In Fig. 1 we show the normalized energy per unit charge E n (X)/n = E n (\) / (Aimr/) for n = 1 
and n = 2 as a function of A G [10~ 5 , 10 5 ]. In the BPS limit A = 0, both curves coincide since 
E n (0)/n = 1 for all n (self-dual solutions). For A > both curves increase monotonically, with 
£^2(A)/2 > £i (A) in agreement with the results of Ref. In the limit of infinitely large Higgs 
coupling constant the curves for n = 1 and n = 2 converge to the constants £1(00) = 1.787 @,H 
and i?2(oo)/2 = 2.293, respectively. Thus the two curves coincide only at A = 0, implying the 
existence of only a repulsive phase. 

In the limit A — > 00 the Higgs potential acts like a constraint, forcing the norm of the Higgs field 

= v/-2Tr(02(f)) = ^4>\(r,9) + 4>l{r,9) 
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to be equal to its vacuum expectation value, = rj. However, this is in conflict with the 

boundary conditions at the origin. Thus, in this limit the norm of the Higgs field approaches 77 for 
all r except for r = 0, where it vanishes. 

For n = 1 the function fa vanishes and the Higgs field is parametrised by a single function which 
corresponds to its norm. In the limit A — > 00, in this case the Higgs field itself is completely fixed 
to its vacuum expectation value for all r except for r = Q. For n > 1 on the other hand, fair, 0) is 
non-trivial and only the norm of the Higgs field is fixed to a constant for all r except for r = in the 
limit A — > 00. Thus the Higgs field itself remains non-trivial in this limit. It can be parametrised by 
the phase function 

$(r,0) = arctan(0 2 (r,6»)/<Mr,#)) • 

As in the n = 1 case, the contributions of the Higgs potential to the energy are negligible in this 
limit, and consequently the energy becomes independent of the Higgs coupling constant. 

The gauge field functions do not develop a discontinuity in the limit of infinitely large Higgs 
coupling constant. Indeed, in this limit they can be identified with solutions of the same model 
with modified boundary conditions for the Higgs field at the origin, which allow for Higgs fields 
with constant norm. We have constructed these solutions with |0(f)| = r\ numerically for n = 1 and 
n = 2. Denoting these solutions by Hi, fa, etc., their gauge field functions H correspond to the 
limiting functions of the gauge fields of the multimonopoles, and, for n > 1, their phase function 
$ corresponds to the limiting function of the phase function $ for all r except for r = 0. Their 
energy E n is finite and coincides with the energy -E n (A — * 00) of the multimonopoles. At the origin, 
however, the energy density diverges and the gauge field functions are not analytio^- Also for n > I, 
the phase function is multivalued at the origin, i. e. $(r = 0,6*) is a non-trivial function of 9. 
Expanding the differential equation for the phase function at r — we find 

$(r = 0,61) = 2arctan(tan"((9/2)) - 9 . 

From this unacceptable behaviour at the origin, we conclude that these solutions are unphysical. 

Let us next consider the monopoles for small values of the Higgs coupling constant. Near A = 
the normalized energy per unit charge can be written in the form 

E n i\)/n ~ l + c„A Q . 

We show E n i\)/n — I for n < 3 and small values of A in Fig. I (inlet). The figure suggests that 
the exponent a does not depend on the winding number n, and that it has the value 1/2 and not 
1 as might have been expected naively from the static Hamiltonian. This result is very natural, 
however, since (A/2) 1 / 2 is just the mass of the Higgs field (in dimensionless units). The physical 
interpretation is therefore, that the increase of the energy with increasing Higgs coupling constant 
A (for small A) is proportional to the mass of the Higgs field. We remark, that the increase of the 
energy for small values of A cannot be obtained by substituting the Prasad-Sommerfield solution 
into the energy functional. In that case, due to the power law decay (~ 1/r) of the Higgs field the 
integration of the Higgs potential would diverge. By considering a Taylor expansion of -E n (A)/n — 1 
in terms of A at A = and taking into account that the gauge fields and Higgs field depend on A 
implicitly, we find that the first term is again the divergent integral of the Higgs potential evaluated 
with the Prasad-Sommerfield solution. This implies, that even for arbitrarily small values of A the 
model with finite A cannot be treated as a perturbation of the model with A = 0. 

Let us now compare with the Weinberg-Salam model of the electroweak interactions, where the 
Higgs field is in the fundamental representation of SU(2). The non-perturbative solutions of this 
model are saddlepoints of the energy functional (sphalerons and multisphalerons). Here sphalerons 



This is also true for n = 1 
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represent the top of the energy barrier between adjacent vacua |l5|,|l3|], and, like multimonopoles, 
multisphalerons are also characterised by their winding number n. Numerical analysis shows, that 
for small values of the Higgs coupling constant, the energy per winding number of multisphalerons 
(n > 1) is smaller than the energy of sphalerons (n = 1), whereas for large values of the Higgs 
coupling constant it is larger, with equality occurring roughly at Mh ~ My/ Because of their 
instability, however, one cannot speak of attractive and repulsive phases of sphalerons. 

Finally we remark that, for n > 2, multimonopoles need not be axially symmetric. Indeed, in 
the BPS limit also multimonopole solutions with only discrete symmetries have been found for the 
higher topological charges p7j . It remains an open challenge to construct such solutions also for 
finite Higgs coupling constant. In particular, a numerical analysis would reveal, which of the possible 
shapes of a multimonopole with a given higher topological charge would be energetically favoured. 
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